Asymptotic behavior of the mean square displacement of the Brownian parametric 

oscillator near the singular point 



Tohru Tashiro 

Department of Physics, Ochanomizu University, 2-1-1 Ohtuka, Bunkyo, Tokyo 112-8610, Japan 

(Dated: March 2, 2009) 

A parametric oscillator with damping driven by white noise is studied. The mean square displace- 
ment (MSD) in the long-time limit is derived analytically for the case that the static force vanishes, 
which was not treated in the past work |l| . The formula is asymptotic but is applicable to a general 
periodic function. On the basis of this formula, some periodic functions reducing MSD remarkably 
are proposed. 

1. INTRODUCTION 

We shall deal with the Brownian motion described by the following Langevin equation which is the same one in the 
previous work [l[: 

x(t) + f3x{t) + [w + #(*)] x(t) = f(t) (1) 

where x(t) is a position of the Brownian particle at time t, /3 is a damping constant per unit mass and f(t) is a 
centered Gaussian-white noise with correlation function 

</(t)/(f)> = 2e5(t-t') . (2) 

4>{t) is a periodic function of 7r whose maximum and minimum values are ±1 as defined in Rcf. 1]. Even if the period 
is not 7r, we can obtain Eq.([T]) again by introducing scaled parameters. In the motion of a charged particle in a Paul 
trap with buffer gas near standard temperature and pressure by Arnold et al. [2J , which is accurately modeled by this 
Brownian motion, w and q in Eq.{T]) correspond to a dc and an ac voltage, respectively. 

It is reported in the previous paper [l[ that for a general periodic function (f>(t), three variances in the long-time 
limit under the condition (x(t)) — {v(t)) = H, 

a x 2(t) = (x 2 (t)) , a v 2(t) = (v 2 (t)) and a xv (t) = (x(t)v(t)) , (3) 

can be obtained analytically by using series of q as long as w > holds. Moreover, it is showed that as q increases 
the mean square displacement (MSD) can be suppressed less than that at q = 0, i.e., MSD of a Brownian motion in 
a harmonic potential. We called this phenomenon a classical fluctuation squeezing. 

In the experiment with the Paul trap, a static force produced by a dc voltage balances with gravity, which means 
that we should set w = in Eq.JT]). For this case with w = 0, however, MSD in the long-time limit cannot be 
represented by a series of q in 

oo 

^ = ^- + V a n q n , (4) 

n— 1 

since the zeroth-order term diverges, which means that (w, q) = (0, 0) is a singular point of MSD in the long-time 
limit. We can understand this from another point of view. Equation {T]) with w = q = denotes a free Brownian 
motion. Therefore, MSD in the long-time limit goes to oo because of a well-known relation, o~ x 2(t) oc t. In short, the 
theory about MSD in the previous paper is right only for the case with w > and cannot be used at w = 0. Thus a 
new method is desired in order to solve the system with w = 0. 

In this paper we will investigate the asymptotic behavior of MSD in the long-time limit near the singular point 
w = q = for general periodic functions, which is not taken up in the previous work. Of course, investigations for 
this case have been done for the Mathieu's case of 4>{t) = 2sin2£ @, but there is no study dealing with the system 
with general periodic functions as far as we know. In fact, we will show that MSD for other periodic functions can be 
smaller than that for the Mathieu's case. This result will reveal a significance of employing other periodic functions 
in experiments using the Paul trap with a buffer gas. 

2. PERIODICITY OF VARIANCES IN THE LONG-TIME LIMIT 

Here, we shall examine the periodicity of variances in the long-time limit which is a very important property when 
investigating the asymptotic behavior of MSD. Of course, the periodicity is clear in the previous paper [1|, but the 
proof is valid as long as w > holds. The proof on this section remains true even if w < 0. 
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Three variances fulfill the following differential equation 



j t X(t) = -{W + q <£(t)]X(t) + K , (5) 



with 

' cr x 2(t) \ I -2 
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(6) 



C7 x0 (t) J \ w -1 /3 

(See appendix [XJ) For simplicity, we represent the coefficient matrix by A(t): 

A(t) = -[W + q$(t)} . (7) 

It is well known that a solution of Eq. ([5]) can be obtained by using the Green's function of the homogeneous equation 
like this: 



X(t) = G{t,t )X{t ) + f dtiG(t,ti)K 



(8) 



where 

n t n t 



pt pt pt\ 

G(t,t ) = l+ dt 1 A(t 1 )+ dh dt 2 A{t 1 )A(t 2 ) + ■ ■ ■ (t > t ) (9) 

J to J to Jto 

in which 1 is the 3x3 unit matrix. 

One can easily confirm that G has the following properties, 

G(t + 7r,t' + 7r) = G(t, t'){t>t') , (10) 
G(t, t")G(t", t 1 ) = G(t, t') (t > t" > t') (11) 

and 

G(t',t) = G(t,t'y 1 (t>f) . (12) 
Therefore, for a natural number N, we obtain 

G(t, t - Ntt) = G(t, t - n)G{t -n,t — 2tt) X • • • X G(t — (N — 1)tt, t - Ntt) 

= G(t,t-7T) N 

= {G(t-TT,0)G(Tr,0)G{t-TT,0)- 1 } N 

= G(t - tt, 0)G(tt, 0) N G(t - tt, 0)- 1 . (13) 

Now, we examine variances in the long-time limit which are obtained by setting t = t — Ntt in Eq.® and then 
letting N goes to oo. If the absolute values of all eigenvalues of G(ir, 0) are lower than 1 which is equivalent to the 
homogeneous solution being stable in the long-time limit, the first term of Eq.® vanishes in such a limit because of 
Eg . (fTB")) . We can find that the second term is also finite only if the eigenvalues fulfill the same condition as follows: 

pt N ~1 pt — n-n 

\ dt x G{t,t x )K = dhG(t,ti)K 

Jt-N-K n=Q 7t-(n+l)ir 

N-l pt 

= X] / dt 1 G{t,t 1 -mr)K 
n =o J*-* 

N-l f t 

= ^2G(t,t-Tr) n dtiG(t,t!)K 

n=0 

— ► {l-G(t,t-7r)} _1 f dhG(t,ti)K (N — ►oo). (14) 

Jt-TT 



3 



In order that this sum of the matrix series converges, it is necessary that the absolute values of all eigenvalues of 
G(t, t — 7r) must be lower than 1. By setting N = 1 in Eq. (fT3")) . one can see that the eigenvalues of G(t, t — tt) and 
G(w, 0) are equal. If the absolute values of all eigenvalues of G(ir,0) are lower than 1, this necessary condition is 
satisfied. 

Then, variances in the long-time limit are given by 



X(t) = {1 -G(M-tt)} 1 / dtiG(t,ti)K 

Jt-TT 



(15) 



By setting 1 1— > t + tt on the above equation one can obtain 

ct+n 



/>t + 7T ft 

X(t+ir) = {l-G{t + ir,t)y 1 dhGit+Tr^^K = {1-G(t,t-Tr)y 1 dtiG(t+n, tx+ir)K = X(t) , (16) 

Jt Jt~ir 

which means that variances have a periodicity of tt. 



3. ASYMPTOTIC BEHAVIOR OF MSD 



At the previous section, we showed three variances in the long-time limit by use of G like Eq. (|15|) . This repre- 
sentation is useful for demonstrating the periodicity of them, but is not practical when investigating the asymptotic 
behavior of MSD. Thus, we shall derive another representation. 

By vanishing o v i (t) and a xv (t) from Eq.([5]) with w — 0, we can obtain 



d 3 cr x 2{t) ^ 3f3 f^(t) 



+ 2{2#(t)+/? 2 } 



dt 3 r dt 2 " l-^^-' 1 <~ > dt 
Since a x 2 (t) in the long-time limit is a periodic function of n, we expand it in a Fourier series 

oo 

a x 2(t)= ]T d n e i2nt (d n =d- n *) . 

71— — OO 

We also expand (f)(t) in a Fourier series 

oo 
rn — — oo 

By putting Eq. fTSj) and Eq. (fT5[) into Eq. fTTj) and arranging it, we can obtain 



E 



4in(/3 + in) (J3 + Yln)d n + Aq {i(2n - m) + 0} 



e i2nt = 4e 



This equation can be expressed as 

A n d n + Aq^2 {K 2n - m ) + P} Cmdn-m = 4e<S„,o (n = 0, ±1, ±2, • • ■) 



where 



A„ = 4in(j3 + vn){ft + i2n) 



and Si j is the Kronecker symbol. 

Let us represent Eq. (|2ip by a matrix formulation 



{A+qB)d = K 



where 



-1 o 



A 



A_i 

Ai 



4(/3-i2)c 4(/3-i)c_! 4/3c_ 2 
4(/3 - i)ci 4/3c 403 + i)c_i 
4/3c 2 4(/3 + i)ci 4(/3 + i2)c 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 
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d = 



I : \ 

d-i 
d 



-l 

o and K = 
i 



( : \ 





V ; J 



-l 
o 
l 



(25) 



We number rows and columns of these matrices and vectors as shown at the sides of them. 
For simplicity, we shall introduce a matrix "D (g) given by 



Hence, d can be solved like this: 



T»{q) = A+qB . 



d = V(q)~ 1 K 



(26) 
(27) 



Furthermore we denote the determinant and the adjugate matrix of T>(q) by A(g) and *D(q), respectively. 
Here, let us estimate o x i (t) in the long-time limit by averaging over the period as past studies [3, 0, Ej| : 

i r r 

<j x 2 = - / dta x z(t) = d (q) = 4e T>(q) 



By using the determinant and the adjugate matrix, the above equation becomes 

4eP , ( 9 ) 

A(q) 



(28) 



(29) 



When q goes to 0, A(q) becomes since all the elements of the zeroth row and the zeroth column of X>(0) are 0, 
which causes the divergence of MSD in the long-time limit. Of course, in order to reach the conclusion it is necessary 
that X>o,0 (0) 7^ 0, but this becomes obvious later. 

We assume that the order of the singular point is n, so that we can expand a^J around q = as follows: 



ct- n a_„+i 

&x 2 = —z- 4 — r ^ h a o + a iQ H 

q n q n 1 



(30) 



It is clear that the coefficient of q n can be obtained by 



4e lim 



:^o,o(g) 



'9^0 A(q) 

If the order is 1, by use of L'Hospital's rule on the above equation a_i becomes 

P ,o(0) io, 
a.-i = 4e — r— — = 4e- 



A(0) 



A(0) 



(31) 



(32) 



in which the over-dot means the derivative with respect to q. 

Note that A(0) is not if cq / 0. This can be understood in the following ways: As is well known, 



A(q) = Tr ±>(q)±>(q) 



Tr 



T>{q)B 



which leads to 



A(0) = Tr [X>(0)B] = Tr \ab] = ^i„, m S m , n 



(33) 



(34) 



where A. is the adjugate matrix of A~ From Eq. ([2^)) . one can confirm that all the elements of the zeroth row and the 
zeroth column of A. are 0. Therefore, the adjugate matrix has only one nonzero element -4o,o[ = 2->o,o(0)], i.e., 



Ao.oS n nSm.O 



(35) 
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Substituting this into Eq. ([34|) . we obtain 



A(0) = Ao,oB Qfi = 4/3c A,o 



From the above equation, we can get 



(3c 



(36) 



(37) 



If Co > 0, <5vT is positive for small q. On the other hand, if Co < 0, 0^2" becomes negative, which cannot be understood 
physically. This contradiction implies that the area where w = and q is small belongs to an unstable region, which 
is made up of parameters which make this system unstable in the long-time limit. Thus, a x 2 cannot be defined. As 
an example, we have shown stable and unstable regions of the system where (f>(t) is a square wave in the past paper 
0. FromFig.2(c) in Ref.Q corresponding to the case with cq < 0, one can see that an unstable region includes a line 
w = within small q. 

If Co = 0, a_i diverges to 00. Thus, the order of the singular point is more than 2. If we assume that the order is 
2, a~2 can be derived as 



tt-2 



A(0) 



(38) 



In this derivation, we have used LTIospital's rule. By differentiating Eq. (|33j) with respect to q and setting q = 0, one 
can obtain 



A(0) = Tr 2>(0)B 



(39) 



The elements of ~D(0) can be calculated as follows: 



£>n,m(0) = < 



■AofiBo,rr, 

A m 
A„ 



■5 n ,o {m ^ 0) 



S m ,o (n ^ 0) 



(40) 



(See appendix IBP It does not matter that £>o,o(0) is not determined by the above equation, since £>o,o = 0- Then, 
equation (|39|) with the above results becomes 

A(0) = ^£ n , m (0)£ m ,„ 



= ^o,m(0)S m ,o + ^«,o(0)B ,r> 

I / \~* Ba^ n B m fl \ ^ Bnfll3o,n \ „ V - ^ ^0,m£>m.O 

- -A),o 2^ — a + 2^ — a — = ~ 2Ao '° 2^ 



A, 



(41) 



Moreover, from Eq. (|24p it is easily confirmed that the elements of the zeroth row and the zeroth column of 3 are 

Bo,m = 4(/3 + im)c_ m and B mfi = 4(/3 + im)c m . (42) 
By inserting these representations into Eq. (|4"Tj) . A(0) becomes 

4(/3 + im)|c m | 2 , ^ P\cJ 2 



A(0) = -2io,0 }Z 



\m{[3 + i2m) 



16io,o 



-± (3 2 + Am 2 ' 



which yields 



Ot-2 



E 



2/3|c m | 2 
m=l /3 2 +4m 2 



(43) 



(44) 
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This means that the order with the case cq = is 2. 
We can reach the following results: 

e 1 



<J X 2 



8 2 J C ™I~ a 



£ 1 1 

(3c q 



^ + 0(q- 1 ) (c = 0) 
0(q°) (c > 0) 



(45) 



Note that the order of the singular point depends on whether Co is or not. 

For the Mathieu's case of <j>(t) = cos2i where |c m | 2 = <5| m |,i/4, the equation ([45]) becomes 

£ 2(/3 2 +4) 

CT z 2 = a 2 r U W J ■ 

8 q- 



(46) 



This result corresponds to the past works0, Q by setting q i — > 2g. 

By these results, we can know the asymptotic behavior of MSD in the long-time limit with a general periodic 
function. These expressions including only one term are useful for finding a periodic function which reduces MSD. 
By comparing the coefficient of q or q 2 , we can find such a periodic function more smoothly than calculating MSD 
haphazardly. 



4. PERIODIC FUNCTIONS REDUCING MSD 



Let us find a periodic function which makes MSD in the long-time limit smaller than that for the Mathieu's case 
of 4>(t) = cos2i in the following two cases: cq = and Co > 0. We set its unit e/8 when showing MSD. 

Here, we review how w, q and <p(t) are determined which is explained in Ref.[l[. w expresses the middle value 
between the maximum and the minimum of an arbitrary periodic function with period tt. q denotes the interval 
between the maximum and w (between the minimum and w). Then, this periodic function is expressed as w + q<fi(t). 
Therefore, the period and the amplitude of <j>(t) must be ir and 1, respectively. Needless to say, cos2i satisfies these 
conditions. We will compare MSD for the Mathieu's case with that for a periodic function satisfying them. 



4.1. periodic function with Co = 



We propose a square wave 



m = { \l {n 7. l) \ t L*j :i )?r ' (» = 0, ±1, ±2, • ■ •) (47) 



(48) 



— 1 if (n — o)7r < t < rwr , 
as a periodic function whose Co is equal to 0. The square of the absolute value of coefficients is 

2{1 -(-!)«} 



k,l 2 = 



With this |c„| 2 , the coefficient of q 2 is derived as 

1 7T/3 3 



v oo 2|c m |2 7r/3-4tanh(7r/3/4) 



(49) 



This coefficient of q~ 2 is smaller than that of the Mathieu's case. We plot two coefficients in Fig[TJ from which one 
can see that the coefficient of a square wave is about half that of the Mathieu's case. Therefore, we can presume that 
MSD for a square wave must be smaller than that for the Mathieu's case. 

In Fig(2l a^2 at (3 = 1 for a square wave and the Mathieu's case are shown. The circles are numerical solutions for 
the Mathieu's case. The solid curve is a solution for a square wave, which we can derive analytically Q. The two 
dashed lines mean asymptotic solutions, Eq. (l45p . for each case. As one can see, both asymptotic solutions correspond 
to the numerical and the analytical solution at a region where q is small, which guarantees that the result in the 
previous section is valid. 
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•s. V, 



-2 -1.5 -1 -0.5 0.5 
logg 

FIG. 2: cT^z at /3 = 1 as a function of q. The circles denote the Mathieu's case (cos2i) derived by solving Eq.([SJ) numerically. 
The solid curve means the analytical solution for a square wave. The two dashed lines are asymptotic solutions [Eq. (|45|l ] for 
each case. 



As we presumed, MSD for a square wave is smaller. However, the two minimum values seem to be close from this 
figure. From numerical computations, it turns out that the minimum for a square wave is 5.15 at q = 1.50 and that 
for the Mathieu's case is 5.25 at q = 1.90. Therefore, the minimum is also smaller. In Fig|3j we show the minimum 
value of a x 2, which we denote by er* 2 and q* which minimizes <j x 2 at several values of f3. The same symbol refers to 
the same /3. From the upper left to the lower right on this figure, /3 increases from 0.5 to 5.0 in steps of 0.5. For this 
range of fj, the minimum values for a square wave are always lower than those for the Mathieu's case. By setting 
logg w 0.4, however, cf^i" for the Mathieu's case can be lower than that for a square wave as seen in FiglJJ which arises 
from the fact that the range of q which makes this system stable for a square wave is narrower than for the Mathieu's 
case. 



4.2. periodic function with Co > 

We define a periodic function with Co > in the first interval from to n as follows: First, we form the positive 
part of this periodic function at < t < T\ from ip(t/T\) where Ti < ir and t^(t) is a non- negative function defined 
on < t < 1 whose maximum value is 1. Next, we form the negative part at Tx < t < 7T from —ip((t — Ti)/T 2 ) with 
X2 = 7T — T\. Then we join these two parts together, so that we can get the periodic function in the interval like 
this: ip{t/Ti) {1 - Q(t - Ti)} - ip((t - T 1 )/T 2 )Q{t - Ti) (0 < t < it) in which 6(f) is the Heaviside's step function. 
Moreover, the 7r-periodic function fulfilling the conditions mentioned at the beginning of this section is completed by 
connecting this segment repeatedly. 
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0.8 
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0.3 



cos It 
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square wave • t # r 



0.2 



0.4 0.6 
log q* 



(3 = 5.0 

0.8 



FIG. 3: CT* 2 which represents a minimum value of o^J vs q* which minimizes 0^2". Gray and black symbols represent values for 
the Mathieu's case (cos 2t) and a square wave, reactively. The same symbol refers to the same (3. As the symbol changes from 
the upper left to the lower right on this figure, /3 gets larger from 0.5 to 5.0 in steps of 0.5. 



Using the periodic function defined in this way, we derive 



c = - / 1 dty{t/T x ) 



TV 



dtip{t/T 2 ) 



Ti-T 2 



dtip{t) 



2Ti 



dtip(t) 



(50) 



Here, we introduce the ratio of the time intervals of the two parts, 7 
function. Then, cq becomes 



Co 



1-7 
1 + 7 



dtip(t) 



T2 /T\ , denoting the asymmetry of this periodic 



(51) 



because Ti = j^— • Since Cq must be positive, it is necessary that < 7 < 1. From this equation, we can find that as 
7 decreases with fixing ip(t), which indicates that the asymmetry becomes extreme, Co becomes larger, which means 
that MSD becomes reduced because of Eq. ([45|) . Therefore, we can reach a conclusion that the asymmetry suppresses 
the fluctuation of position for the case w = as well as the case w > reported in Ref. [l[ . 

As an example of such a periodic function we show an asymmetric cosine function £(t) (see Fig]?]) which can be 
obtained by setting ip(t) = sin7rf and inserting the negative part into the positive part divided equally. Of course, 
this procedure does not affect the above discussion. Thus, we can get 



2 1 



C() 



7 



7T 1 



Moreover, we find 



(1 



-i2nT, 



11 )2i (1 + 



)72i 



4n 2 T- L 2 



7T 2 - 4n 2 7 2 T 2 



(n = ±l,±2,--.) . 



(52) 



(53) 



Figure [5] shows 0^2 at (3 = 1 as a function of q. The triangles and the boxes represent numerical solutions for an 
asymmetric cosine at 7 = 0.1 and 7 = 0.5, respectively. The circles represent numerical solutions for the Mathieu's 
case. Dashed lines represent asymptotic solutions described by Eq. J45|) for each case which coincide with numerical 
solutions at small q. Since the orders of the singular point are different between the two cases, these lines are crossing. 

For small q, the difference between MSD for an asymmetric cosine and for the Mathieu's case is extreme. MSD 
at 7 = 0.1 is smaller than that at 7 = 0.5, which can be easily understood as follows: From Eq. (|52|) . cq decreases 
monotonically as 7 increases. The coefficient of g _1 is proportional to 1/co, and so MSD increases when 7 increases. 
For small q, the ratio of MSD at 7 = 0.1 to at 7 = 0.5 is about 40.7%. The minimum value of MSD for the Mathieu's 
case is larger than that for 7 = 0.1 and 7 = 0.5. Even if 7 gets bigger, the large-and-small relation does not change: 
from numerical calculations, it is clarified that although 7 varies up to 0.9 the minimum for an asymmetric cosine 
does not become more than that for the Mathieu's case at the range of j3 from 0.1 to 5.0. 




FIG. 4: An asymmetric cosine function The time intervals where £ is positive (7\) and £ is negative (T2) are not equal. 
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FIG. 5: op" as a function of q at /3 = 1. The triangles and the boxes represent values derived from Eq.([5]l numerically for an 
asymmetric cosine at 7 = 0.1 and 7 = 0.5, respectively. The circles represent the Mathieu's case (cos2t). Dashed lines are 
asymptotic solutions described by Eq.(gSJ for each case. 



5. CONCLUDING REMARKS 



In this paper, we considered analytically the Brownian parametric oscillator described by Eq.([T]) especially with 
setting w — which is not treated in [l[ and the following points are made clear: I. It is proved that the period of the 
variances becomes the same as that of 4>{t) in the long-time limit. II. asymptotic behaviors are derived as Eq. (|4"5)) . 
which are the main results of this paper and are applicable to general periodic functions. III. A square wave and an 
asymmetric cosine are proposed as periodic functions with cq = and cq > 0, respectively, which make MSD smaller 
than that for the Mathieu's case. 

One of main results, ct^T = e/(/3coq) + O(q ), can be obtained more easily from another point of view: If q is very 
small, the motion can be separated into a "fast" part whose period is it and a "slow" part which hardly changes 
during the period [6j]. By taking the time average of the total motion, the fast part vanishes. The effective potential 
which induces the slow part is written as 0] 

2 00 I |2 \ 

^o + ^E^Rif, (54) 
2 * — ' m j 

m=l / 

where x(t) denotes the slow part of the motion. Hence this slow part is a Brownian motion in the harmonic potential 
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and its MSD in the long-time limit can be calculated like 



whose first term (proportional to QL" 1 ) is compatible with our result. On the other hand, equation (I55j) with cq = 
does not correspond to our result [9(. If we assume that (3 is small, however, they become equal. Therefore, it is the 
achievement of this paper that the asymptotic solution for cq — can be acquired precisely. 

Our results indicate that MSD can be reduced more remarkably than ever in experiments where this Brownian 
motion can be observed, e.g., experiments using the Paul trap with a buffer gas 0. The author hopes that this result 
will be beneficial in such experiments. 



6. ACKNOWLEDGMENT 



The author would like to thank Dr. Akio Morita for the extensive discussions. 



APPENDIX A: DIFFERENTIAL EQUATION GOVERNING THREE VARIANCES 



The derivatives of three variances shown in Eq.f3|) are 

{£)) \ a v2 (t) + (x(t)v(t)) 



d 
dt 



a x 2(t))\ ( 2a xv {t) 

(a At)) = 2(v(t)v(t)) I . (Al) 



From Eq.([T]), we can obtain 

(v(t)v(t)) = -(3(v(t) 2 ) ~[w + #(*)] (x(t)v(t)) + (v(t)f(t)) (A2) 

and 

(x(t)v(t)) = -0{x(t)v(t)) -{w + q<f>(t)} (x(t) 2 ) + (x(t)f(t)) . (A3) 

{v(t)f(t)) and (x(t)f(t)) are derived by use of Novikov's theorem [4( which states that a functional of Gaussian noise 
g[f{t)\ satisfies 

(<?[/(*)]/(*)> = f <#mw)) ( s -j0y) 

where Sg[f(t))/5f(t') indicates the functional derivative of g[f(t)} with respect to /(£'). 
Because of 

x(t) = x(0) + f dt'v(t') (A4) 



and 



we get 



which yield 



v(t)=v{0)+ f dt'{-/3x(t')-[w + q<f>(t')]x(t')+f(t')} , (A5) 



Sx(t') 6v(t>) 

I7W) = ;0and W) = l - (A6) 

(x(t)f(t)) = and (v(t)f(t)) = e . (A7) 



Therefore equation (|A1[) becomes 



{x{tf) \ / 2 

£ | (v{tf) \ = \ -2(3 -2[w + q <l>(t)] i i (v(t) 2 ) i - | 2e | . (A8) 

(x(t)v(t)) \-[w + qcj>(t)} 1 -d 



dt 



(v(tf) - 


(° 


h 2e 


x(t)v(t)) j 


\ 
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APPENDIX B: DERIVATION OF V n , m (0) 

As is well known, T>(q) and T>(q) satisfy the following relation 

T>(q)V(q) = V(q)T>(q) = A(q)E (Bl) 
where E is the infinite unit matrix. By differentiating this relation with respect to q and setting q = 0, one can obtain 

T>(0)A=-AB (B2) 

and 

AD(0) = BA . (B3) 

In this derivation, we have used A(0) = 0. 

The elements of the nth row and the mth column of both hand sides of Eq. (|B2j) are 

i>nMAm = A m P„, m (0) (B4) 

i 

and 

- Y^AnjBLm = - ^ Aofi8 n $8iflBi, m — -AofiB , m 6 n fi . (B5) 
i i 

Keeping Ao = in mind, we can derive 

Km(0) = ~ A °f°' m S n , (m ^ 0) . (B6) 
Similarly, from Eq. (|B3[) we can also derive 

KnM = - Ao f n >° S m . (n^O). (B7) 
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